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non-lin&Ttrc dc plasmas cm configuration de queue rnajinctiquo, comme dans lc cas <le la 
mugnclusphere terrestre, Finalcment nous utilisons de precedent resultats pour presenter unc 
methods auln-coherente d" analyse dc la xtnbilitt non lincairc du plasma cl do hi transition 
coire lcs LOinporteiPcnf linc-airc et non I infinite pour des systcmcs prcs tic rcqullibrc. La 
methods se base sgr I'analysw de la density d'energic poteiUiclle, utilisant des rdsultais 
ric emplacement du plasma provenant d'un mode'le Iin6airc pour caiculcr tss energies mix 
dcuxiemc ci tmisicme ordres. Nous dfimontrons la mc'thode dans un cxcmple dc champ lihrc 
avee des lignes de champ dipolaim Altrees. Dans ecu exemple, nous voyon* cluircmem la 
transition emre lea instabilitcs lin£aires et non Uncaires. 

|Traduit par In Redae(ion| 

1. Introduction 

The problem of nonlinear plasma instabilities is of a great importance in many areas of space 
and laboratory plasma physics. These instabilities have many common characteristics even though the 
conligurations of the plasmas might be quite different I J | s and are characterised by [21 

■ sudden onset, 

* buret of energy, momentum, and particles across magnetic surfaces, 

• avalanche character (localized instability can disrupt the whole system ]3 1). 

For example, in space physics, observations show that during the onset of the inlensiJication of a magneiic 
substorm a huge amount of energy Is suddenly released on time scales of the order of minutes (Alfvcnsc 
time scale) |4|. These time scales suggest that the processes involved are highly nonlinear (explosive). 
There were several previous attempts to create a working model of explosive instabilities in magnetized 
fluids 13, 51. Approaches based on perturbation theory in the framework of magnctohydrodynamics 
(MHO) typically lead to extremely complicated expressions, 

A very practical way to address the question of nonlinear stabilities appears to be a variational 
approach |6-8|. The advantage of a variational approach h that one is dealing with a scalar quantity 
(energy), which makes this formalism suitable in any magnetic-field geometry. 

Wc present an alternative mathematical approach to the Lagrangian-Hamiltonian method, employ- 
ing differential geometry [9-1 1], using Tensor notation combined with general transformation rules for 
the change of coordinates. Further, we use the generalized Stokes theorem to derive equilibrium con- 
ditions and wave equations from the expansion of the l^igrangian. This approach is very effective in 
eliminating redundant terms, and in minimizing the complexity of the algebra. In addition, the results 
are independent on the topology of the magnetic field. 

Our derivation follows the general method fur the expansion of the Lagrangian of ideal MHO 
outlined in rcf. K. Treating the plasma displacement as a coordinate change, wc derive expressions for 
the magnetic field, plasma pressure, and density up to the third order in plasma displacement. Then we 
express the ideal MHD Lagrangian u p lo the third order in plasma displacement using these expressions. 
In the third-order term, wc correct an error in the contribution of the magnetic field in rcf. 8. We then 
derive equilibrium conditions and general wave equations for plasmas without ambient convection 
from the first- and second-order Lagrangian term, respectively. This linear equation makes possible 
calculations of plasma eigenmodes in complicated topologies of the magnetic field, since it is derived 
with no restriction on the topology of the field. Also, we present several suggestions on the use of this 
formalism in the investigation of nonlinear plasma instabilities. Wc also show the consistency of our 
results with the established works on two special configurations of the plasma. 

Further, we apply certain restrictions to the wave equation and derive the equations for MHD 
field-line resonances [1*2] and linear ballooning modes [13], and we demonstrate that these results are 
compatible with published results. These examples allow us to illustrate the strengths of our approach 

^ 2004 NRC Canada 
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and make a disconnection between our results and those using vector. formalins. However, we show 
that out method is self-consLstcnt and much more transparent than the vector approaches. 

S^^P^ • ™ del tbat in^igatcs whether certain configuration, of magnetized plasma, 
can lead to an onset of explosive behavior. Our method is based on a comparison of the second- and 
third-order potential energies for the plasma displacement obtained as a solution of the linear equation. 
The assumption behind this method is that the system we are investigating is not far from equilibrium, 
and that initially linear behavior, whether stable or unstable, is dominant. During the temporal evolution, 
however the system can undergo a transition from linear to nonlinear behavior. We demonstrate the 
use of this method with an example of a magnetic-field topology that is stretched from the dipolar 
configuration similar to the stretching in the magnetotail during the growth phase of the magnetosphcric 
substorm r 14]. Since the stretching of the field lines provides a source of free energy in the system, due 
to enhanced currents f I K evolution of an instability might be expected in this system. The model wc will 
present here does not attempt to describe the full dynamics of the i notabilities, since these dynamics are 
beyond the scope of ideal MHD, which assumes sufficiently slow processes and large spatial scales. 

Our work has four objectives. The first is to develop a formalism thai will simplify what arc now 
very complicated calculations using variational methods in plasma physics. Tn this regard the methods 
of differential geometry have been very effective in other areas of physics where complex topologies 
might arise (see, for example refs. 10 and 15). Geometrical methods have also been used in connection 
with plasma physics 116, 17]. Wc then address two further objectives, the correction of the results in 
ref. 8, and an outline of methods to solve some problems in plasmas with complicated magnetic-field 
topologies. Last, we present an example of the nonlinear stability of a cold plasma in a stretched field 
topology, 

2. Lagrangian In ideal MHD 

Although ideal MHD gives one of the simplest approximations of plasma dynamics, jt still describes 
some very interesting physical features of the plasma. This relative simplicity gives us an ideal tool for 
developing new methods for studying the behavior of plasmas. In ideal MHD, the Lagrangian has the 
simple form (the notation is given in the Appendix) 



(1) 



This Lagrangian needs constraints to yield nontrivial behavior. The constraints we use are mass con- 
servation f v pdV - const, magnetic flux conservation f s B . &S = const and the adiabatic condition 
pp-Y = const. 

We treat plasma perturbations as a coordinate change 

i m X + 0 (9 

where the displacement £ is a function of tbe original coordinates *, and time f. Also, we assume that 
the displacement £ and its gradients (f j) are small enough so that we can use them as the ordering 
quantities for a perturbation treatment of the plasma. Then the whole pertu rbation calculus can be done in 
a geometric framework, without introducing any new physical constraints on the system. The geometrical 
approach uses a component formalism, combined with the use of definitions for the transformation of 
the volume and surface (instead of deriving them through additional physical assumptions [8], and 
uses geometrical definitions of determinants and inverse matrices during the derivation process. This 
formalism simplifies the algebra significantly, and the results are in a compact form that is suitable 
for further calculations. It also has the advantage that the results do not depend on the coordinate 
system. Consequently, the results hold for any geometry, and do not limit the physics beyond limitations 
introduced by use of ideal MHD. 

02OO4NRC Canada 
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2.1. Transformation of kinetic energy 

The perturbed velocity has the same form v = dx/di as the unperturbed velocity, with the difference 
that the plasma displacement is also an explicit function of x and. therefore, total different .atoms 
necessary Taking into account mass conservation, and the fact that dV = JdV. which implies that 
p = p/J. where J is the Jacobian of the coordinate change (2). the perturbed kinetic energy has the 
form 

In the absence of ambient convection the first and the third terms in brackets vanish and the kinetic 
energy is simply 

ic^jdvlpm) 2 (4) 

This will be the form for the kinetic energy used throughout the rest of the work. 

2.2. Transformation of thermal and magnetic energy 

Since mass conservation combined with the rule for the transformation of a volume element yields 
p = p/J % the application of the adiabatic condition, connecting density and pressure yields the trans- 
formation relation for the pressure in the form p = pJ~ Y . ^ 

Similarly, since the surface element dS transforms as &St 2= A/dS/, as will be shown in the next 
section together with the derivation of the transformation matrix A, magnetic flux conservation yields 
B' = (A-^YjB J . So the potential energy terms are 

for the thermal component, and 

E M -/ 4vlj 8lJ (A"') V (A-')V (6) 

where gtj is a metric tensor, for the magnetic energy. Therefore, combining (4) with (5) and (6), the 
perturbed Lagrangian in the absence of convection is 

£ = / a v(l, ( a, s) =-^l-i^(.-«)V(A-')V) m 

So far our results arc consistent with those derived in ref. 8. However, our method of derivation of 
the Jacobian and the matrix A differs significantly. The method we will present in the next section is self 
consistent, purely geometrical, and does not require any additional physical assumptions, for instance, 
restrictions on the magnetic field. The use of the component formalism and the use of generalized 
Kronecker deltas and Lcvi-Civita symbols ensures a very compact form for the results, making further 
calculations much simpler. Also, the authors of ref. 8 made the error of assuming that the inverse matrix 
A -1 is second order. Although this ta true for the matrix A, jt is not true for the matrix A -1 . Due to the 
omission of these terms, their expression for the third-order perturbation in the Lagrangian cannot be 
generally correct In the special case they discussed, with the trivial third-order term, these additional 
terms vanish as well. At the end of this section, we present a more detailed analysis of the terms that 
ref. 8 omitted due to the missing third-order term in the expansion of A" 1 . Another important point 
is that a major advantage of the formalism wc use occurs in situations with special symmetries. Many 
terms disappear due to their symmetry properties, significantly reducing the complexity of calculations. 

0 2004 NRCCnnada 



Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 

PAGE 22/34 1 RCVD AT 2/28/2005 1 1 :01 :21 AM [Eastern Standard Time] 1 SVR:USPT0€ FXRF-1/6 1 DNIS:8729306 * CSID: 1 DURATION (mm-ss):0948 



597 

Dobia3 and Samson 

2.3. Transformation of volume and surface 

Starting with the definition of theJacobian J = det(^) and u«in e the expression for the determinant 
of an n x n matrix A, det A{ « S{n)A ••■A [9, 1 1], the expression 

j = (4 + ei) (57 + «3) w + «) * 1 + «i + j + «w*J*'«i (8) 

for the J acobian is obtained. , _c 

The surface element is defined as dfi - « y * dV dx* . Using this definition for the perturbed surface, 

nodng that d# = {4 + $j k ) dx*. and using the fact that the product of symmetric and antisymmetric 

tensors is 7.ero, the expression 

d* = Mjr (4 + d) (< + *i) **° dx" * Aj iSt W> 
for the transformation of surface elements is obtained. This also defines the transformation matrix A 
used in the previous section. However, for the calculation of the magnetic energy density, the inverse 
matrix >% -t is needed. This inverse matrix has a form similar to the original matrix A. and 

(a- 1 )'. = 3s%„ «■ + er' w +er' (,0) 

The inverse matrix {J~X = («f + $.*)'' is by definition expressed as (J-^ s= ty)W£ n J?J?. 
Expanded to the third order it is 

Substituting (I I) into (10) aad expanding the result to the desired order one obtains the inverse matrix 
A~ x with arbitrary accuracy. Expansion up to the third order yields 

(a-^). = 6) + ^ - - + \ &> l {M + *&) 

Equivalent results for tbe transformation of the magnetic field could be obtained in this special case 
by using the Lundquist identity, which can be written in the form B ! = B* (&) + V tJ ) I J 1 1 8]. Instead, 
we derived this relation directly from the magnetic flux conservation to demonstrate this method in its 
most general sense, with the direct use of the constraint to Lagrangian (1). 

2-4. Perturbed Lagrangian up to the third order in plasma displacement 

How it is possible to substitute expressions for tbe Jacobian (8) and (12) into expressions for 
the thermal and magnetic energy (5) and (6). 

Thermal energy 

Substituting expression (8^ ( 1 +J) ] ~ y = 

l + (j-y)j-(I/2)y(l->/)y 2 -r-(l/6)y(l-y 2 )/ 3 in(5)one obtains tbe expression for the thermal 
component of the potential energy up to the third order 

*«/ d^[«+(i-v)(^+?M+^ 1 ) 

-!> (#J + itftftf ) + ^ CJ3) 

©7004 NRC Canada 
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2,4.2, Magnetic energy 

After substituting expressions (8) and (12) into expression (6) one obtains for the magnetic energy 

+ so*'*' («j (i - *5 + -i- - «i*U?) < l4 > 

The magnetic energy density consists of two components, the first one is magnetic pressure and the 
second one takes into account curvature. In atypical vector formalism for deriving (13) and (14). there 
is often a very large number of terms [8], and it is often not clear which terms are zero and which terms 
cancel. Our use of component notation eliminates this problem, giving far fewer terms, and equations 
that are much easier to use in the analysis of stability problems. 

2.5. Expansion of the Lagrangian 

Using expressions (4), (13), and (14) it is possible to expand Lagrangian (7) with respect to the 
plasma displacement as 

The terms £° to £ 3 are 

P — --B 2 (16) 

y - I 2 

C°> = (p + ±* 2 ) # - gijB'B'^, (17) 



i 

2" 



+• Jay e^'^fjC - (*j*S*5 +fctt5«i - (19) 

The sr,eroth-ordcr term is simply the unperturbed Lagrangian of ideal MUD. Applying the variational 
condition / dr / d VCf^ 0 to the first- and second-order terms, tl is possible to derive the equilibrium 
condition and momentum equation, respectively. In the next section, we will discuss the linear wave 
equation in more detail, and we wilJ show consistency between our results and other established works. 
The third-order term contains information about possible explosive behavior 18]. This is discussed in 
more details later. 
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2 6 Effect or the third-order contribution in transformation of magnetic field 
" As was mentioned earlier, the authors of ref, * omitted the third^order term in the derivation of 
the matrix A' 1 that is used in the transformation of the magnetic field. They used this matrix as an 
intermediate step in tbeir derivation of the matrix A % which is indeed of the second order in p1™»^ 
placement. However, A" 1 also contains higher order expansion terms and the^e terms heeome important 
in the transformation of the magnetic licld ^ 

Since a direct comparison between the expression from the Ugrangian obtained in ref, 5 with 
expression (J 5) is very difficult to perform due to the complexity of their result (containing o large 
number of redundant terms), we decided instead to trace the difference that omission and (or) inclusion 
of the third-order term in the expansion of A" 1 causes in the final expansion of the magnetic energy. 

A" 1 appears in the magnetic energy u> the expression gy {A"^ B k BK The magnetic 

pressure term in ref. 8 would be (3/2) (j^ft^tS). whereas our term is (1/2) (B^%fe%) . This 
introduces an asymmetry between thermal and magnetic pressure terms in Pfirseh and Sudan's results. 
This will affect any system with nonvanishing V • 

All the other terms affect the curvature part of the magnetic energy. In fact, the only curvature 
term that remains with no third-order term in the A -1 matrix is simply (1/2) fa B k B ! S'fytty and 

the whole expression g u B* B } + titj# " * omitted. This means that any system 

with the plasma displacement changing in any direction other than along the field lines will not be 
described correctly. Again, in the situations with special symmetry that the authors of ref. 8 use to verify 
their result, the tbird-order term is identically zero due to the symmetry properties. In these special 
instances their results and our results agree. 

3. Variational calculus 

Using a variational approach in the second (17) and the third term (18) of the Lagrajigiau (15), it 
is pos$ible to derive equations governing the dynamics of the plasma in a linear approximation, The 
condition /d V>C {1) = 0 yields the equilibrium condition in the form 

-(p + 5* 2 ), + *'^ s j-° (20) 

This is the standard equation containing the sum of thermal and magnetic pressure and pressure due to 
field-line curvature. 

The condition f d* fdV£( 2) = 0 yields the linearized momentum equation for plasmas in the form 

+ tfBi (^»,^(*'a' - W)) - (**0 fl 3 " 7 * - * PbJ ) 

+ (%«-w/)(^B^(e , *'-*y)) wi) 

Obviously (21) is a wave equation describing linear wave propagation in a plasma with arbitrary (and 
possibly very complex) magnetic-field topology, for which idcAl MHO is a valid approximation. There 
are no further restrictions on the direction of propagation nor on the configuration of the ambient 
quantities. The only restriction on the validity of this equation is that the wave amplimde be small 
enough for the linear approximation to be valid. However, by applying restrictions on this equation, it is 
possible to derive dispersion relations for waves in specific plasma configurations. Since the methods of 

fi> 2004 NRC Canada 
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differential geometry are not commonly used in MHD, in the following section WO give ^J*™**" 
SdStrion^ 

that equation (21) yields result consistent with established results using .standard vector methods, but 
in a more transparent way. 

Xi. Field-Hne resonances (FLRs) 

In the case of a rectilinear magnetic field the natural choice for coordinates is a Cartesian coordi- 
nate system. Then the metric tensor gij = diag(l, 1.1). and vectors and covectors are identical. The 
equilibrium condition (20) takes the simple form (P + B*H\i = «■ IWog this in equation (21), and 
assuming a harmonic time dependence, the momentum equation has the form 

-(BjB*)9p(f J B'-pB') (22) 

We now choose axes such that the z axis is along the magnetic field lines, and the variation of the 
ambient quantities is along the x axis. The magnetic field is B = (0, 0. B) and all derivatives of P, B, 
and p except 3 X vanish. Also assuming that the displacement is harmonic along the y and z directions 

t{p9„tix) = |W expi(A,y + k.z) 0 s * 
equation (22) can be reduced to 

P (" 2 - *2 vf) «, = -a, (/ ? (x)3.,b) < 24 > 
where f is defined as 

y? = *I, v? = ^. v»-vj + vj. * = 

P P 

This result is consistent with the expression obtained in ref. 12. Equation (24) yields two turning points 
and two romances. In the case of a cold plasma, expression (25) reduces to 

This equation has only one turning point at a> 2 k 2 V* and one resonance at a> 2 = k\ V% t Occurrence 
of this resonance is due to coupling between the compressional Alfvcn wave and shear AJfv6n wave of 
the same frequency. Since there is no curvature, ballooning modes cannot evolve and thus this mode 
cannot be explosively unstable. Tn the last part of thus paper, we will show energy trnasfer in this mode 
as part of the demonstration of the stability investigation. 

Though wc presented a simple model of FLRs in a Cartesian topology, the strength of our method 
will be more evident in complex magnetic-field topologies. One example of more complex topologies 
occurs in the region of stretched field lines in the near-Earth raagnetotail [14) during the substorm 
growth phase. Though these modes play an extremely important role in the plasma dynamics of the 
magoetospherc, we know of no complete solutions (analytic or numerical) of the equations for the 
cigenmodes. The linear wave equation (21 ) is written in a form that is suitable for such a calculation. 
AH that needs to be defined is the metric tensor of the topology we are working in. If this metric is 
connected to the magnetic field, then all the information we need is contained in this tensor. 

©2004 NRC Canada 
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3.2, Linear ballooning mode 

Wc consider a second example with a curved magnetic field and an ortbonormal coordinate system 

defined by 

„ B (27) 

i - . (29) 
<t> = cxn 

Here again, the metric is gij = diag(l. ». D> and vectors and coveetors arc identical. Introducing the 
definitions 

K h = /I.VIn B (31) 

Kc =«-(|.V)? (32) 
the equilibrium, condition (20) can be written in components as 

8^-0 < 33 > 

3^=0 ( 35 > 
Assuming no azitnuthal dependence (fy = 0) of the ambient field, the problem reduces to two 
dimensions. Equation (21) can be split into components in the parallel and radial directions. Assuming 
that the ambient magnetic field changes slowly along the field lines, and consequently assuming a 
harmonic dependence of the displacement along the field lines (i.e., 3j| ify), and in time (i.e., 
d t -ia)), the parallel component of equation (21) is simplified to 

(36) 

Subsequently, assuming very small-scale sizes in the azimuthal direction (k+ -y oo), which is equivalent 
to averaging the total perturbation pressure I\ = &p + B to zero [131, which in the terms of 
equation (21 ) is 

Pr = yPd^ + S'SiP - B,d p -$**') = 0 W 

and using (36), the equation 

fvj + V 2 _ 9/t J = _2K c Vfa . (38) 

connecting the congressional and shear Alfven wave, is obtained. 

Tn the absence of the perturbation pressure, the radial component of equation (21) is simply 

|y - v}h\ - avfoc* + *)]«, - 2r e V* - ^ J Brf' (39) 

This equation describes the generation of the shear Al tven wave by the compressional wave. Equa- 
tions (38) and (39) describe a linear ballooning mode, and are consistent with results obtained in rcf . 1 3. 

O 2004 NRC Canada 
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3.3. Cold plasma approximation 

As the last example, we shall derive an instability criterion in ihc cold plasma approbation. This 
approximation Climes plasma pressure P = 0. The.v using an orlhonormal coordinate system with 
coordinate along the magnetic held, the wave equation (21) reduces to the form 

. t n m 

pfcb - JL (mm* + m**) + ft ^> 

where A r are the Lame coefficients (square roots of the diagonal elements of the metric tensor). Since 
the parallel displacement £i is constant, wc can generally assume it is equal to /.cro and the problem 
reduces to two equations. Assuming a harmonic temporal and spatial dependence of the displacement, 
with infinite azimulhal wavelength, and assuming thar the ambient quantities vary much more slowly 
than the plasma displacement we obtain the dispersion relation for the plasma 

Hcrc v 2 _ and /? 2 a«e Lame coeflicicnts with coordinates *| andx 2 .and*i and k 2 are wave 

vectors in directions xi and jc 2 « respectively. liquation (43) yields two solutions tor (o 1 



wf _ v ?f!I-5l (44) 



- 2 = ^§ (45) 



This solution allows for instability if 

4 < 4 (46) 

In the magnctosphcrc-like configurations, the wavelength along the Held line is often much larger than 
the radial" scale size, these systems can easily yield unstable behavior, More detailed discussion of 
instabilities in force-free field can be found in rcf. L 

4. Explosive instabilities 

4.1, Outline of the method 

The third-order term in the l-agrangian (15) can, under certain circumstances* provide us with the 
information about the nonlinear behavior of the plasma, without the need to understand the nonlinear 
dynamics of the system. There are some instances (some of the configurations with special symmetries 
that have the third-order energy identically zero) when this method will not work, We will not be 
considering these special configurations, instead we assume that the expansion (15) contains all the 
terms. 

The method we outlined is based on the observation that equations like 

x<xx* W) 
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wh icb correspond to * mechanical system with the potential energy V * x\ have singula solutions of 
the form 

(t - t 0 r 

xv.Av.fort, we ^sumc that systems with dominant third-order potential energy (or, in other words, a 

Systcrn the third^der potential energy starts to dominate, we arc looking at the onset ot exp os.ve 
behavior. This idea was outlined in ref. 8. however, the complexity of the expressions makes it very 
difficult to apply this idea to specific problems. The presence of higher order terms can possibly lead 
to the saturation of the instability [19]. This is, however, beyond the scope of our discussion here. Our 
goal is to identify the onset of the instability, not to resolve the full dynamics of the instability. 

4.2. Instability investigation 

As we have outlined, analysis of Lhc energy can give us stability properties of the plasma in the 
nonlinear regime without having to solve nonlinear equations. The whole process can be dmded into 

three steps. . , , . . 

First, it is necessary to define the equilibrium or quasi -equilibrium configuration of the plasma that 
we want to investigate. This will define our ambient pressure, magnetic field, and plasma density P, B s 
and p. The magnetic field. B< then defines the geometry of the problem. 

For the second step, we solve the linear wave equation (21 ) to obtain the class of allowed displace- 
ments for a gryen configuration. Obviously, it would he useless to investigate all displacements, since 
most of them would be nonphysical anyway. We assume that the solution of the linearized equation 
of motion is a good approximation of the real displacement, since we do not attempt to describe the 
nonlinear stage, we only want to predict its beginning. This approximation is justified a posteriori by the 
dominance of the second-order energy at the initial stage of the model. At this stage it is also possible 
to derive the linear growth rate and thus the linear stability criterion. 

In the third step, we evaluate the second- and third-order potential energies and compare them. When 
the third-order term start* to dominate, the system has reached the threshold of explosively unstable 
behavior. If the second-order term dominates, linear behavior prevails. This linear behavior is already 
resolved in step two. During the stage when the second- and third-order terms are comparable, we are 
looking at a transition stage with the nonlinear dynamics being important in the system, and, thus, using 
the linear equation for the description of the plasma is not correct. However, if this transition stage is 
short, and then evolves back into either a linear stage or the third-order term becomes dominant, we can 
still predict the behavior of the system. 

As was mentioned earlier, the fourth-order term can contribute to the saturation of the instability. 
Since we are not identifying the mechanism of the saturation, it is not necessary to use the fourth-order 
terra in the present analysis. However, to understand the full dynamics of the instability, the fourth-order 
term becomes very important and needs to be accounted for [19]. 

5. Explosive instabilities in stretched magnetic-field topology 

To demonstrate the investigation of a possible explosive instability, we decided to use a configuration 
of the magnetic field that would be similar to the Earth's raagnetosphcre in the slightly excited state 
on the night side just prior to the explosive release of energy in substorrn intensifications (Fig. 1 ) [ 14] . 
Although the plasma in the near-Earth plasmasheet is high 0, we used a cold-plasma approximation 
since an initial equilibrium (J x B = 0) is easier to find and wc shall be using this example only for 
illustrative purposes. The topology of a force-free field is more complicated than the hot plasma case 
in Earth's magnetosphcrc. In our configuration, we stretch the field so thai the projection of a field 
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Fie 1 Comparison between a dipolar topology and the Stretched field line topology we wed. Field line 
It 10?Ki to 12Re at thTequmor. Number, denote corresponding field lines. Th.Sju projection 
of the stretched field lines on die meridional plane. In the force-free configuration these field bnes hove an 
ayjmuthnl component. 

5 



h 0 




Fig. 2. Dependence of oj 2 on the coordinate (along field lines). The frequency is normalized to the order 
of J. Clearly, in the region close to the equatorial plane the system is unstable. 
2 

f 1 

JS 



%.2 ™4U5 -O-OS 0 0.15 0.2 

X 1 



line on the meridional plane has a configuration similar to those shown in rcf. 14. Since the condition 
j x B = 0 requires the current to be parallel to the magnetic field, the magnetic-field line wil) not 
be in meridional plane, but must have an azimuthal component. The cold-plasma equilibrium can be 
found by the use a connection between the magnetic- field geometry and coordinates (fij B/B), In 
other words, the magnetic field can be derived from the coordinates as B ^ M/g\\. Therefore, for a 
cold plasma, it is possible to define the coordinates lirst. and then simply calculate the magnetic field. 
From the conservation of magnetic flux (V ■ B = 0) it follows that the relation #n = gizgx must hold 
for the coordinates used. Then we solve (40)-(42). Solutions of these equations arc then substituted 
into the expressions for the second- and third^ordcr potential energy. The solution of (40) is trivial, thus 
reducing the problem to two dimensions. 

Assuming that the ambient magnetic field changes on a larger spatial scale than the plasma displace- 
ment, the dispersion relation for this system is approximated by (44). with the Alfven velocity defined 
as v i = jpip ^ an d k 2 are the components of the wave vector. Figure 2 shows an example of the 
dispersion relation (44) in the case of equal wave length in both the radial direction and along field 
lines. This dispersion relation allows instability in regions close to the equatorial plane. However, due 
to the scaling of the problem, it is more reasonable to expect the radial wavelength to be shorter than 
the wavelength along a field line. In such a case the area of the instability can be greatly enhanced. 

For the numerical solution of (41) and (42). wc used the method outlined in rcf. 20- We used 
normalized units with the unit of length equal to 1 Earth radius. Time was normalized so that Alfv6n 
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ir„ r>vn am i« of the system in the equatorial plane (top), and the second- and third-order energy 

Tmt !h, nondimeJonal unit, The linear tartfi, ^ 
appears to saturate, then it restart? again 

0.15 



and eventually the energy shows onset of explosive behavior. 




1000 



velocity in the equatorial plane would be 0.01, The magnetic moment was 0,05. The grid has 256 
(along field lines) x 128 <io radial direction) grid points. The time step dr «- 0-00001 and we ran 
for 0 < t < 1000 which allowed sufficient time for possible interesting features to develop. The 
ambient magnetic field and plasma density were defined as B = M//t t , where we put M = 0.05 and 
p t= w/ sin 8 9, with a> = 10~ 5 . We had to use such a small value for the magnitude of the ambient 
magnetic field due to numerical stability problems. The azimuthal wave number was k* = 2. The 
ionospheric boundary was put at 4 (Earth radii) to avoid problems with significant gradients at the 
boundary. The initial plasma displacement and velocity were 

5 f . 0 = 0, i=/i.v,0 < 49 > 
,1,0 = 0, J-p,* (50) 

/ (r-r 0 ) 2 \ ( (9 - ft) 2 \ 

where we put A = 0.0005, r 0 = 13, d = 0.2. 6b = *A &n<* * - 0.0J5- This models the initial 
perturbation of the initial equilibrium as a very small impulse in the radial direction. We varied the 
topology of the field from dipolar to various levels of stretching. As expected, test showed the dipolar 
configuration to be stable. 

Two examples of the \ nstabili ty are shown in Pigs. 3 and 4. The top of each figure shows the dynamics 
of the system in the equatorial plane. In both cases linear instabilities arc observed. The bottom of the 
figures shows the behavior of the energy densities* using the expansion of the Lagrangian ( 1 5) calcul ated 
using the results of the linear equations (40H42). We can observe that at the time corresponding to the 
saturation of the linear instabilUy the thitd-order term becomes dominant. This indicates a break down 
of the l inear approximation and the onset of an explosive phase. 
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Yin 4. Dvmtmlcs of the system in the equatorial ptaM (top), "nd the second- and third-order ™<f 
critics (bottom). Time i* in nondlmariiaal units. After short linear -stage the system becomes oxpkwvcly 
unstable 
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6. Discussion 

6J- Instabilities in coW plasmas 

A question can arise concerning the nature Ol'Lhc instabilities in coW plasmas Ufcc the model con- 
sidered here. At first we need to look at the source of possible fvee energy that can drive the instability. 
Such a source of energy is the parallel currents ,/|| that produce the stretching (and twisting to give an 
a/.Smuthnl component) of the magnetic-field lines in this model of cold plasmas. Tn Tact, it is this energy 
connected to the stretching that is the primary source for the instabilities even in case of high-beta 
plasmas as was demonstrated in rcf. 2 1 . 

Tn the case of cold plasmas, the equilibrium condition is simply J x B = 0. This expression can 
be alternatively written for ideal MHD as VB 2 /2 = B VB by using the fact that J V x B. This 
condition can be compared with the equilibrium condition for a warm plasma, V{P + B 2 /2) - BVB. 
In this case, perturbations of the equilibrium between the total pressure (combined plasma and magnetic 
pressure) and the curvature can generate a ballooning instability in the plasma. An excellent discussion 
of Raylcy-Taylor instabilities in space plasmas is given in ret', 22. It is reasonable to assume from the 
similarity of the equilibrium conditions between cold and hot plasma*, where in the former case the 
magnetic pressure takes the role of the total pressure, that the nature of the insUibility in a cold plasma 
will be effectively the same as in case of a hot plasma with the magnetic pressure gradient acting against 
the field line curvature to generate ballooning. That way it is possible to have an instability even in the 
absence of a thermal pressure gradient in a plasma. 

6.2. Topology of the instability 

In a force-free magnetic field, the equilibrium condition j*B = 0 requires the current to be parallel 
to the magnetic held. This leads to the existence of the azirautbal component of the magnetic field to 
allow for stretching of the field lines io the meridional plane. Thus, the magnetic-field line is not entirely 
in one plane, but is twisted across the meridional planes. 
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If we allow for a small displacement of the field line &D Earth-ward, a current is generated in the 
azimuthal direction. This current combined with the ambient j t generates a magnebc-Held component 
g ho azimuthal component of the ambient field. Reduction of the a^muthal component ofthe 
fie d causes in turn the generation of current feeding the Earth-ward displacement ofthe held lines. Thi* 
way the instability evolves, and eventually bring* the field back to a dipolar state. The energy that was 
stored in enhanced currents is transformed into the kinetic energy of the enhanced flows in a plasma. 

7. Conclusion 

This paper has addressed several objectives. First, we have shown that Ac use of geometrical 
methods such as component formalism, transformation relations, and the formulas for determinants and 
inverse matrices can simplify the description of magnetised plasmas. Compared with a traditional vector 
approach, the derivations arc more transparent, fully self-consistent, and the results are in a compact 
form due to the easy reduction of redundant terms. Also, in systems with special symmetries, tr^se 
symmetries are relatively easy to incorporate into the equations. 

We have applied this method to calculate the expansion of the Lagrangian within ideal MHD, but the 
same method can be applied to oonideal MHD models by adding additional terms to the unperturbed 
Lagrangian, or alternatively, the method can also be used for other approaches to plasmas with a 
known Lagrangian. We also present suggestions for the use of this formalism in plasma problems that 
are complicated by nontrivial topologies of the magnetic field, such as calculation of eigenmodes in 
asymmetric topologies, or investigations of nonlinear instabilities in curved magnetic-field topologies. 

Treating plasma perturbations as a coordinate change, a method common for elasticity theory (231, 
we expanded the Lagrangian of ideal MHD up to the third order in plasma displacement. The first two 
orders of the expanded Lagrangian yield equilibrium conditions and the general linear wave equation 
for that model. We have shown that for ideal MHD this equation is consistent with the results obtained 
by the direct linearization of the MHD equations using a traditional vector approach. We have also 
derived dispersion relations and a linear instability criterion for a cold plasma in a general curvilinear 
field. The third-order term provides information about nonlinear plasma stability. In systems close to 
equilibrium, we can use the results ofthe linear model to calculate the second- and third -order energy 
densities. Their relative size provides insight into the initial stages of nonlinear dynamics. Specifically, 
when the third-order energy prevails, the system has evolved into an explosively unstable stage. 

Last, we demonstrated this method in the investigation of a cold plasma in a curved, stretched 
magnetic-field topology where the results reveal the possibility of explosive behavior that follows an 
initial linear instability. In our study with a cold plasma model, the energy source for the instability is 
contained in the stretching ofthe field lines beyond the dipolar configuration and the current enhancement 
connected with this stretching. This two-stage instability is consistent with the results of nonlinear MHD 
modeling [24]. This qualitative agreement suggests that this method is indeed an effective tool for the 
identification of the onset of explosive instability. 
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Notation 

B ambient magnetic field 

fj ambient pressure 

v velocity 

p plasma density 

Y isotropic coefficient 

Ay metric tensor 

A determinant of the metric tensor f> {j 

t; f k fully anti-symmetric unit tensor 

S/mir generalized Kroneeker delta (product of 5s anti symmetrized in all indices ) 

a u -j] antisymmetrization in i, j indices 

a, ; = 'M 
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